Achieving the threshold regime with an over-screened Josephson junction 
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We demonstrate that by utilizing an over-screened Josephson junction as a noise detector it is 
possible to achieve the threshold regime, whereby the tails of the fluctuating current distribution are 
measured. This situation is realized by placing the Josephson junction and mesoscopic conductor in 
an external circuit with very low impedance. In the underdamped limit, over-screening the junction 
inhibits the energy diffusion in the junction, effectively creating a tunable activation barrier to the 
dissipative state. As a result, the activation rate is qualitatively different from the Arrhenius form. 



In the course of scientific progress, it is desirable to 
use recent advances of physical understanding to develop 
new experimental techniques that will in turn give impe- 
tus for further advances. This description is apt for the 
current state of electron counting statistics in mesoscopic 
physics, and the aim of the current Letter. Recent ex- 
periments in electron counting statistics have measured 
the asymmetry of the current distribution [1, 2], and the 
effect of the measurement circuit [1, 3]. Subsequent re- 
searchers have taken several independent approaches to 
exploring this physics. If the system is transferring elec- 
trons on sufficiently long time scales, it is possible to 
directly count individual electrons [4], look at various 
higher current cumulants [5] , and even examine the con- 
ceptually new area of conditional counting statistics [6] . 
It is also possible to focus on the frequency dependence of 
the cumulants [7, 8]. A further exciting possibility is to 
measure the tails of the current fluctuation distribution 
via a threshold detector [9, 10]. 

A natural threshold detector is a Josephson junction 
(JJ): by measuring the rate of switching out of the 
metastable supercurrent state into the running dissipa- 
tive state, information about the statistical properties of 
the noise driving the system may be extracted. This sys- 
tem has been the subject of recent experiments [11-13]. 
However, it was noticed that a single JJ typically works 
near the Gaussian point [9] (a regime that has been well- 
studied in the past, see e.g. [14, 15]), so much so that 
the third current cumulant makes only a small correc- 
tion to the escape rate [16-18]. This is because the slow 
semiclassical dynamics of the JJ averages out the Marko- 
vian noise source and becomes effectively Gaussian under 
typical conditions. The purpose of this Letter is to show 
that despite this difficulty, the threshold regime is realiz- 
able for a single underdamped J J. The threshold regime 
goes beyond the third cumulant and realizes the full po- 
tential of the Josephson detector, where the escape to 
the running state of the JJ is driven by the tails of the 
distribution, rather than by relatively small deviations 
from the average. 

Circuit effects are known to be important for electron 
counting statistics. The measurement circuit can lead to 
cascade corrections [19, 20] to higher current cumulants, 
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FIG. 1: (Color online) The negative logarithm of the escape 
rate is plotted versus dimensionless barrier height AU /Ec for 
various values of the threshold parameter V ■ The solid lines 
correpond to Poissonian noise, while the dashed line corre- 
sponds to a QPC with transparency T = 0.06. Inset: Simpli- 
fied electrical circuit for the Josephson junction (J J, marked 
with an X) threshold detector of the noisy mesoscopic system. 

masking the system contribution. A similar effect has 
been predicted by the authors to occur in JJ detectors 
[17]. The central idea of this Letter is to have a very 
small load resistance, resulting in an over-screened junc- 
tion, whilst still maintaining the underdamped state (see 
inset of Fig. 1). The advantage is that not only cascade 
corrections are suppressed but also the relative contribu- 
tion of the third cumulant is enhanced compared to the 
Gaussian contribution [17], assisting in the measurement 
of the third cumulant in recent experiments [12, 13]. We 
will demonstrate that this is the case not just for the 
third cumulant, but that all higher cumulants are also 
enhanced, leading to non-perturbative activation by rare 
events. 

Josephson detection circuit. — The circuit in the inset 
of Fig. 1 shows the essential part of the detector com- 
prised of the J J with Josephson energy Ej , and the ca- 
pacitor, C. The fluctuations in the circuit originate from 
the combination of the current II through the macro- 
scopic load resistor and by the mesoscopic system current 
Is, which is to be measured. According to Kirchhoff's 
law, the total current ^5 + II is equal to the sum of the 
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Josephson current 1$ = (Ej/&o)sm(j) where $o = H/2e, 
and the displacement current Iq = CV. This leads to 
the equation of motion for the superconducting phase </>, 



C^l4> + Ejan<f> = 9 (Is+lL), 



(1) 



where we used the relation V — &ocf>. 

We consider an ohmic system and load resistor, so that 
(I s ) =J S - G S V and (I L ) = J L - G L V, where G s , G L 
are the system and load conductances, and the constant 
currents Js = GsVs, Jl — GlVl are tunable parameters 
that will be shown to control the activation threshold. 
Equations (1) can be rewritten as a set of Hamiltonian- 
Langevin equations for the phase variable 4> and canon- 
ically conjugated momentum p = $oQ (where Q = CV 
is the total charge on the capacitor), 



= p/m, p = -dU/d(j> + $ 5I, 



(2) 



with "mass" m — &qC, and where SI = Is — Js + Il — Jl 
is the dissipative part of the system and load current. 
Equations (2) describe the motion of a "particle" in the 
tilted periodic potential 



U(<f>) = -Ej[coB<f> + J<t>], 



(3) 



stimulated by the dissipative part of the system and load 
current, where J — &o(Js + Jl)/ Ej is the dimensionless 
total current bias. Dissipation leads to relaxation of the 
JJ into one of its metastable supercurrent states, where 
the phase is localized in one of the potential wells so that 
(V) = 0. In the dissipative state, the phase drifts along 
the bias which generates a non-zero voltage drop V. 

Weak damping threshold limit. — Here the phase oscil- 
lates with the plasma frequency uj p \ = — J 2 ) 1 ^ 4 , 
where fij = yjEj/vn. The energy relaxes slowly with 
rate (Gs + Gl)/C < u> p i to the local potential mini- 
mum. We further assume the separation of time scales, 
max{eV~s,T} > fkv p \, so that the noise source Is is 
Markovian. According to our previous results [17], the 
escape rate F (predicted to exponential accuracy) is 



logT = -$ 1 J XdE, (H{\4>)) E = o. 



(4) 



Here, the function H(z), which we refer to as a Hamil- 
tonian, generates the cumulants (irreducible moments) 
of SI in Eq. (2), so that d™H\ z=0 = ((SI n )). The func- 
tion X(E) is the escape trajectory in the extended en- 
ergy phase-space, or the "instanton line" . The notation 
(. . .)e = (l/7p) § dt(. . .) denotes time-averaging over a 
physical trajectory in the (q,p) phase-space at a certain 
energy E, with T p (E) being the period of the quasi- 
periodic motion at that energy. See Ref. [17] for a de- 
tailed derivation. 

We note that after replacing the node voltage V with 
$o0j the first two coefficients in the expansion of Tt are (i) 
(SI) = -(G s + G L )$ O and (ii) ((<5/ 2 )) = 2G L T + ((If)). 



We now focus on the limit of small load resistance, 
Gl S> Gs so circuit back action can be neglected [17]. 
In this case, the load controls coefficient (i). However 
we observe that for a large system voltage, the system 
shot noise will dominate the load resistor noise in co- 
efficient (ii) if eV s > G L T/G S - The circuit load, be- 
ing a macroscopic resistor, has no higher current cumu- 
lants, so the Hamiltonian takes the form (H(\</>))e = 
— Gl$ o A(0 2 ) e+Js(F(\<P))e, where F generates the gen- 
eralized Fano factors of the system, 



F( Z )-f-^ 



(5) 



The instanton equation (4) may be brought to dimen- 
sionless form 



E C X (F(A0))j 



((A^) 2 )i 



eG L 



(6) 



where Eq = e 2 /2C is the capacitor's charging energy, 
and we define the threshold parameter V, which is the 
ratio of the J J energy relaxation rate Gl/C to the meso- 
scopic system's excitation rate Js/e. If relaxation is 
weak, ?< 1, then Eqs. (5) and (6) may be truncated at 
the first (Gaussian) term, giving the escape rate (4) as 



logic 



AU 



(7) 



where T c s — ((/|))/2Gl is an effective noise temperature, 
so the Arrhenius form is recovered in agreement with [17]. 
However, if the relaxation rate dominates the excitation 
rate so ? > 1, then we enter the threshold regime where 
the solution to (6) is non-perturbative in the cumulant 
expansion. 

The physical reason for this is that while the Gaus- 
sian contribution to the system noise cannot compensate 
for the fast energy relaxation, rare current events (where 
many electrons are sequentially transmitted in a short 
amount of time) may be able to excite the JJ strongly 
enough to overcome damping. The power of these cur- 
rent kicks is proportional to the velocity <fi, and therefore 
creates positive feedback for rare-event activation in the 
weak damping regime, because the velocity continues to 
grows as the particle gradually ascends the potential well 
to the escape point. 

By introducing the J J quality factor Q = uj p \C/Gl > 
1 and the separation of time scales parameter 1Z = 
eVs/huj p \ > 1 we reformulate the condition for the 
threshold regime as V = {1 / QK) (e 2 / 4hG s ) > 1, there- 
fore Gs < e 2 /h. This implies that the system is a tun- 
nel junction, which is known to create Poissonian noise 
(however, see the discussion below). The circuit back ac- 
tion can be neglected if eVs > eV = e<f>o0 ~ \JEqE. 
We will show below that E <~ Ec at the threshold, so 
eVs/Ec = QU(hG L /e 2 ) > 1. Together with V > 1 this 
gives the overscreening condition Gl > Gs- 
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Poissonian noise. — In order to illustrate the thresh- 
old behavior, we consider a simple harmonic potential 
U{4>) = (l/2)ma>pj(/) 2 with a sharp cut-off at <f> = </> . 

In this simplified model <j)(t) — y / 2E/Ej cos (oj p \t) . The 
averages in Eq. (6) over the periodic orbit at constant 
energy may now be done exactly. The Poissonian noise 
generator is F(z) — e z — z — 1. Using the integral rep- 
resentation of the modified Bessel function of order 0, 
Io(z) = (l/27r) Jq 277 <i#exp(— zsin0), we find Eq. (6) sim- 
plifies to 



h(z) - 1 





(8) 



This equation is numerically solved for A and integrated 
over the energy E. The result for the logarithm of the 
escape rate is plotted in Fig. 1 as a function of AU for 
different values of V . 

For Pel, the Bessel function may be expanded as 
I (z) = 1 + z 2 /4 + ■■-, whicn gives A = 2hV/E c . When 
A is substituted into (4), the result (7) is recovered, be- 
cause for Poissonian noise ((/§)) = GsVs- In the oppo- 
site limit, P>1, the right hand side of Eq. (8) is large, 
and wc utilize the asymptotic form of the Bessel func- 
tion, Io(z) — > e z J 1 \phvz~ for z large. Thus the variable z 
must be only a logarithmically growing solution of Eq. 
(8): z = L H (E) = \og{2Vyj2irE/E c ), up to the double- 
logarithmic correction. This gives A = Ln(E)^m/2E, 
so the escape rate (4) of activation due to Poissonian 
noise may be approximated by 



logTp = - 



AU ( l2irAU\ 



The dominant square-root behavior is clearly seen in Fig. 
1, and is a sign of the break-down of Gaussian noise ac- 
tivation. 

For the escape to be achievable, we require that AU is 
not much larger than Ec- More rigorously, for strong bias 
the action scales as ^ AU / E c - (1 - J) 3/4 ^/Ej/E C - 
The total number of states in the quantum well, N tot ~ 
AU /fiwpi, also scales down as N tot <~ (1— J) 5 ^ \f~Ej /E c , 
but must be large in the semiclassical limit. Therefore, 
we estimate — logTp ~ iVtot log V j\J\ — J s t . In Refs. 
[12, 13] this number varies from 10 to 17. This indicates 
that it may be hard but nevertheless feasible to observe 
the threshold behavior. 

Anharmonic correction. — Here we show that a real- 
istic potential U (4>) leads to only a small anharmonic 
correction to the escape rate plotted in Fig. 1. In the 
case of a general potential, the instanton line Eq. (6) 
may be found by noting that for Poissonian noise, the 
integral in (F(X4>))e m the threshold regime is expo- 
nentially dominated by the largest value of <f>. En- 
ergy conservation, (m/2)cb 2 + U(4>) — E, indicates that 
this will be near the bottom of the potential, where 



the potential is approximately harmonic with frequency 
lo p \. Therefore, the previous Bessel function result 
will hold to excellent approximation. The other av- 
erage, (4> 2 }e 1 is generalized by noting that by using 
dt = d<j)/\(f)\ and conservation of energy, the equation 
d E [T p (E)(4> 2 } E ] = T p (E)/m holds. This equation may 
be integrated to find (<p 2 ) E = 2irhN(E)/mT p (E), where 
N(E) = (l/27rft) dE'T p {E') is the number of quan- 
tum states in the cavity below the energy E. 

Combining results for (F)e and {c/) 2 )e and going to 
the asymptotic threshold regime, we conclude that the 
instanton solution takes the previously found form A = 
L(E)y / m/2E with the logarithm replaced by L(E) = 
Lh(E) - log[A(E)], where the function 



A(E) 



T P E 
2nhN 



Ed E \ogN 



(10) 



characterizes the anharmonicity. It takes the value 
A = 1 at the bottom of the potential well and di- 
verges logarithmically at the top of the barrier as A ~ 
\og[AU/(AU — E)]. Thus this factor may compensate 
the large parameter V in close vicinity of the barrier top, 
AU — E ~ AUe~ v , so the Gaussian noise there domi- 
nates. However, the overall anharmonic contribution to 
the large logarithm L(E) is relatively small and can be 
neglected. 

Stabilization effects. — In what follows we wish to illus- 
trate an important fact: statistics of rare current fluctua- 
tions carries complementary information about a Marko- 
vian process which is not contained in any finite current 
cumulant. As a first example, we consider the shot noise 
from a quantum point contact (QPC), which is known to 
be a binomial process. The noise Fano-factors are gen- 
erated by F(z) = (1/T) log[l + T(e z - 1)] - z, where T 
is the transmission of the QPC. The cumulants are ob- 
tained by the expansion (5) at z = 0. Therefore, in the 
tunneling regime, Tel, the logarithm can be expanded 
to lowest order in T so the noise of the QPC is Poisso- 
nian, F(z) — e z — z — 1, as in the example considered 
above. However, our results indicate that the Poissonian 
process can always provide a strong enough current fluc- 
tuation for the JJ to escapes from the metastable state 
with some small probability. This is not the case for a 
QPC in the tunneling regime, as we demonstrate below. 

Indeed, the rare current events of the QPC are deter- 
mined by the asymptotic of F(z) at large z: 



1 U + logT, z-> +oo, 
f(Z) + 2= Hlo g (l-T), ^-oo. 



(11) 



This result implies that with a small probability T M the 
current acquires its maximum value J max = Js/T = 
e 2 Vs/irh, so that all M = eVsAt/nh electrons arriving at 
the QPC during time interval At are transmitted. Simi- 
larly, with the probability (1 — T) M all the electrons are 
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reflected giving zero current [21]. If AU exceeds some 
critical value, the maximum or minimum current fluctu- 
ation creates insufficient bias for the JJ to escape from 
the supercurrent state [9]. We named this phenomenon 
the Pauli stabilization effect because it originates from 
the Pauli principle of electron occupation [10]. 

To estimate the value of the bias, J s t, and escape rate, 
T st , at the stabilization point, we average F in (11) over 
the energy-conserving trajectories, which yields (F) e = 
(1/T) {AA^/Tp + (1/2) log[T(l - T)}}, where A</> is the 
distance between the two turning points at energy E. 
The denominator of (6) is the same as before, so the 
instanton line is given by 



X s t(E)= aTp 



Acjy-bN 1 



(12) 



with the coefficients a = -(1/2) log[T(l - T)] and b = 
16nVT. Integration of this instanton line leads to sin- 
gular behavior when Acf> = lQirVTN . Evaluated at 
E = AU this equality determines the stabilization point 
Jst below which the rate T vanishes. We find this point 
by assuming the strong bias limit 1 — J <C 1 in (3) giving 
a cubic potential: U/Ej = (1 — J)<p — 3 /6. Skipping a 
number of steps, we present the result: 



n _ 7 ^3/4 _ Sy/^/gJ 



(13) 



Next, we note that at the critical point J = J st the 
integral (4) with A = X s t(E) from Eq. (12) is convergent 
because A st (E) has only has an inverse square-root diver- 
gence at E = AU. Therefore, logr st can be estimated by 
dropping the least singular term, bN, in the denominator 
of A s t (E) and replacing the potential with the harmonic 
approximation U = (l/2)mu)^ l (p 2 . Straightforward eval- 
uation then gives logr st 
strong bias AU 



y/AU/E c log[T(l-T)]. At 
(1 - Jst) 3 ^ 2 Ej, and using the result 
(13) we find the escape rate right before stabilization: 



io g r st 



log[T(l - T)\ 
VT 



(14) 



We estimate (14) in terms of the quality factor Q > 1 
and the separation of time scales parameter 1Z > 1 as 
logT st ~ QR\og[T {1-T)\. Alternatively, using the total 
number of states N tot > 1 in the quantum well, and Eq. 
(13), we find that lo g r st ~ N tot log[T(l-T)]/VT=J7 t , 
which agrees with the previous estimate for Poissonian 
noise. To compare with the Poissonian case, the escape 
rate for a QPC with the transmission T = 0.06 is plotted 
in Fig. 1 in the case of the simple harmonic potential with 
a sharp cut-off. The sharp potential leads to a logarith- 
mic divergence at the stabilization point AU/Ec = 6.86, 
in contrast with the rate discontinuity discussed above 
for a smooth potential. 

Our second example is the telegraph process: The sys- 
tem current switches randomly from the value I\ to the 



value I2 and back with the rate 71 and 72, respectively 
We have found in Ref. [22] that the cumulant generator 
of the telegraph process is given by the formula 



H s (z) =Iz-i+ y/(AIz- A 7 ) 2 /4 + 



7i72, 



(15) 



where Al = 7 2 - h, A7 = 72 - 71, I = (h + h)/2, 
and 7 = (71 + 72V2. In the case of slow switch- 
ing, 71,72 <C ii,/ 2 , the noise becomes super- Poissonian 
with the current distributed between I\ and 1^- The 
sharp cut-off of the distribution function at these val- 
ues results in the stabilization effect. The asymptotic 
form Hs = Iz — 7 + \AIz — A7I/2 at \z\ — » 00 leads 
to the result (12) with the new coefficients a = 7/ AI 
and b — %-kGlEq / 'eAI . Therefore, at the stabilization 
point the results (13) and (14) hold after the replace- 
ment VT -» G L E c /2eAI and log[T(l-T)] -» -2j/AI. 
With the new separation of time scales requirement, 
VJ = 7/Wpi > 1, the action at the stabilization point 
can be estimated as — logr st <~ QTZ' , so the telegraph 
stabilization effect should also be observable. 
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